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We present a measurement of the e'^e KgRi cross section in the energy range 
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' ^/s = 1.04 1.38 GeV. For the energy ^/s > 1.2 GeV the cross section exceeds vector 
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] taken into account. Measured cross section agrees well with previous measurements. 
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I. INTRODUCTION 

The spectroscopy of light-quark vector mesons is still far from completion. It is mostly 
because not all the processes are measured with sufficient accuracy. One of the sources of new 
information on vector meson spectroscopy is a process e"'"e~ KsKl, into which contribute 
isoscalar u;(783), 0(1020), and isovector p(770) resonances as well as their higher mass 
excitations. Accurate measurement of the e~^e~ — >■ KsKl cross section is also important 
because it is a part of the total cross section of the electron-positron annihilation into 
hadrons, which enters into calculations of the hadronic vacuum polarization contribution to 
the anomalous magnetic moment of the muon and running electromagnetic coupling constant 
at Z-boson mass aem{Mz)- 
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At present the e~^e~ — > KsK^ cross section is measured with an accuracy of few per cent 
only in the narrow energy interval around 0(1020) resonance ^jQl- First measurements of 
this reaction at higher energy were conducted with DM1 and OLYA jj] detectors. The 
DM1 measurement at DCI collider covers the energy range ^/s = 1.40 -r- 2.18 GeV with 
a total integrated luminosity of 1.4 pb~^. The OLYA measurements in the energy range 
^/s = 1.06 -ir 1.40 GeV were carried out at VEPP-2M collider with the total integrated 
luminosity of 0.7 pb~^. In both experiments significant excess of the e~^e~ KgK^ cross 
section over vector meson dominance (VMD) model was observed^The latest measurement 
of the e"'"e^ — > KsKl cross section was done by CMD-2 detector ^ in the energy range ^/s 
from 1.05 up to 1.40 GeV. 

In this work the measurement of the e^e~ KsKl cross section in the energy range 
from 1.04 up to 1.38 GeV is reported. The experiment was carried out with SND detector 
at VEPP-2M e^e~ collider. The analysis is based on the total integrated luminosity of 
9.1pb~\ 



II. DETECTOR AND EXPERIMENT 

The SND detector ^1 operated at VEPP-2M e^e~ collider complex from 1995 up to 
2000. It was designed as a universal detector for the studies of the decays of p, cu, 
resonances as well as the processes of e^e~ annihilation into hadrons in the energy range 
v/i = 0.40 ^ 1.40 GeV. 

The main part of the SND detector is a three-layer scintillation electromagnetic calorime- 
ter consisted of 1632 counters with Nal(Tl) crystals. The total thickness of the calorimeter 
for the particles originating from the detector center is 13.4 radiation lengths. The energy 
resolution of the calorimeter for photons is cte/E = 4.2%/ E{GeV). The angular resolu- 
tion is equal to cr,^, cr^ — 1.5°. The solid angle coverage is 90% of 47r. The charged particle 
tracks are measured by ten-layer drift chamber system located inside the calorimeter. 

In this analysis the data collected in the experimental runs of 1997 and 1999 are used. In 
1997 two scans of the energy range from 0.96 up to 1.38 GeV with an energy step of 10 MeV 
and the total integrated luminosity of 6pb^^ were performed. The year 1999 scan was done 
in the energy range from 1.02 to 1.34 GeV with a step of 10 MeV and a total integrated 
luminosity of 3.1 pb~^. In this analysis the points with ^/s > 1.04 GeV were considered. 
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III. DATA ANALYSIS 

The process 

e+e- KsKl (1) 

was studied in the decay mode Ks vr^vr^ 47. The Kl mesons, due to their large decay 
length, which is much larger than the the detector radius, and large nuclear interaction 
length in the Nal(Tl) (~ 0.35 m), do not produce any signal in the detector in a significant 
part of events. Nuclear interaction of the meson or its decay inside the detector produce 
energy deposits in the calorimeter counters, which are interpreted by the event reconstruction 
program as one or more photons. In both cases appearance of charged particles in the 
tracking system is improbable. Thus, for the analysis of the process (0) only the events 
with no charged particle tracks in the drift chamber are considered. The main background 
processes for the process under study are the following: 

e+e- ^ un^ ^ ttOtt^, (2) 



e+e ^ 0(7) ^7(7) ^ 37r°7(7), (3) 

e+e" ^ 07 ^ KsKlJ (4) 

The process dl} is a radiative "return" to the resonance due to emission of photon(s) by 
initial particles. These photons are mostly emitted at small angle with respect to the beam 
direction and thus are not detected. The process Q is a sum of the processes e+e~ rjj 
and e+e~ — ^ 7777 with additional photon emitted by initial particles. Other background 
sources considered in the analysis are beam background and cosmic particles. 
Initial event selection is based on the following criteria: 

- -^7 > 4, where is a number of reconstructed photons; 

- Nc = 0, where A^^^ is a number of reconstructed charged particles; 

- events with a cosmic particle track reconstructed in the calorimeter are rejected. The 
track in the calorimeter is a group of calorimeter crystal hits positioned along a com- 
mon straight line. 
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The last requirement reduces the number of selected events by more than a factor of two, 
almost completely rejecting the background from cosmic particles. 

The events satisfying these criteria are kinematically fitted in the hypothesis of Ks — »• 
TT^TT*^ 47 decay. The kinematic fitting procedure searches for the combination of two 
photon pairs from 7r° decays. The invariant mass of found vr^'s is constrained to the mass 
of Ks meson. There is no constraints on Ks energy. In a multi-photon events for all four- 
photon combinations corresponding are calculated, and the combination with a minimum 

(Xxs^27ro) is chosen. The 

Xxc^27r0 distributions for different selection criteria are shown 
in Figs, n and El For further analysis the events with xi:g^27rO ^ ^5 are selected. The 
following additional selection criteria are applied to these events: 

1. < (i = 1..4), where Q is a "quahty" parameter of a reconstructed photon equal to 
— log L, where L is a likelihood function value of a hypothesis that observed transverse 
energy distribution in the cluster of hit calorimeter crystals corresponds to one isolated 
photon [7J. This parameter provides separation between events with isolated photon 
showers and those with merged showers or clusters from Ki-meson nuclear interactions 
or decays. 

2. 36° < 9i < 144°, where 6i is a polar angle of the photon included in the reconstructed 
Ks meson with respect to the beam direction. This criterion rejects significant part 
of the beam background. 

3. 400 < Mrec < 550 MeV, where M^ec is a recoil mass of the reconstructed Ks meson: 

Mrec = ^J^s{^S-2Eks)+MI^, (5) 

^/s — total energy in the center-of-mass frame, Eks — the energy of the reconstructed 
Ks meson, is a iiT^-meson mass. This requirement suppresses background from 

the process (jH). 

4. x^o^o^ > 60, where x^o^o^ is a of the kinematic fit in the e^e" -k^-k^^ hypothesis. 
This criterion is applied to events with A^^ > 5 to suppress the background from the 
process Q. 

The total number of events satisfying all criteria described above in the full energy range 
is equal to 1998, of which 585 events are in the region ^/s > 1.1 GeV. 
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FIG. 1: The XKg^2nO distribution for the 
events with ^/s > 1.2 GeV, A'^^ < 7, and Etot > 
0.5 • y/s. Dots with error bars — data, hne — 
simulation of the process (0) , shaded histogram 
— simulation for the process 



FIG. 2: The XKg^2nO distribution for the 
events with ^/s > 1.12 GeV and Etot < 0.5 • ^/s. 
Dots with error bars — data, hne — simulation 
for the process shaded histogram — exper- 
imental distribution for the beam background. 



The number of background events of the process e~^e~ un^ was estimated from Monte 
Carlo simulation. The cross section of the process (0) was taken from the paper Q], in which 
was also shown that the simulation of the multi-photon events in SND detector agrees with 
experiment within 5%. We checked the accuracy of the e~^e~ — > ujtt^ background estimation 
using experimental data. The error estimation is based on approximation of the Xks-*2tt0 
distribution (Fig. H]) for events with y/s > 1.2 GeV, < 7 and Etot > 0.5 • y/s, where 
Etot is a total energy deposition in the calorimeter. These requirements effectively suppress 
contributions from all background processes except (j2)). The distribution is approximated by 
a sum of the process under study and 0, obtained by simulation and shown in Fig.^ This 
experimental estimation agrees within statistical accuracy with that obtained by simulation. 
The total number of events of this background process in the full energy range with all 
selection criteria applied is A^'^^^o = 11.3 ± 0.3 ± 2.3. 

The number of background events of the process e~^e~ — 777(7) is also estimated by 
simulation. The events of this process satisfying all selection criteria are dominated by 
radiative return to (p resonance: e^e" — > 07, —>■ rij, in which the photons emitted by 
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initial particles are not detected. The accuracy of the estimation of this process contribution 
is ~ 3%. It is determined by the accuracy of the e^e' —>■ rj'-f cross section measurements 
near the resonance (~ 2%), and the accuracy of the e+e~ rj'-f'-f cross section calculation 
(1%). The total number of the background events of this process is 21.2 ± 0.2 ± 1.1. 

For the estimation of the beam background the x'ks^2tt0 distribution shown in Fig.|21was 
used. The dots with error bars show the experimental x'Kg^2-K0 distribution for the energy 
range ^/s > 1.12 GeV with additional selection criterion Etot < 0.5^/s, which rejects events 
of all background processes except beam background. The shaded histogram shows the 
experimental distribution for the beam background events, for which the selection criteria 
inverse of 1 and 2 were used. The total number of the beam background events is estimated 
as A^b = 30 ± 3 ± 5. 

The contribution of the e+e^ KsK^j events strongly depends on the beam energy. 
For the energy region close to resonance the photon energy is small in comparison with 
the total energy ^/s making this process virtually indistinguishable from e'^e^ KsKl. 
With the increase of collision energy the photon energy grows and the kinematics of the 
processes and becomes more distinct. The recoil mass spectra against reconstructed 
Ks meson are shown in Fig. El for four energy intervals. The peak with a mean value close 
to the K°-meson mass is due to to the reaction e~^e~ — KsKl, the rest corresponds to 
the process e+e~ — > KsK^'y. Good separation between these processes can be achieved at 
energies above 1.2 GeV. At lower energy the processes and Q cannot be separated. To 
solve this problem the approximation of the cross section was carried out with the detection 
efficiency as a function of both ^/s and the energy of the photon emitted by initial particles. 

The distribution of the number of selected events, background from e"''e^ — cuvr" and 
e~^e~ —>■ rj'-f {j) as function of energy with the beam background subtracted is presented in 
Table □ 



IV. DETECTION EFFICIENCY 



The detection efficiency for the process under study was determined using Monte Carlo 
simulation. The simulation takes into account photon emission by initial particles P. Il0|. 
permitting us to take into account the dependence of the detection efficiency e{^/s, z) on 
energy y/s and z = the fraction of energy carried away by the photon emitted by 
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K^-meson recoil mass, MeV/c 



FIG. 3: The recoil mass spectra against reconstructed Ks meson in four energy intervals: a) 
^= 1.04- 1.05 GeV, h) y/s = 1.06- 1.09 GeV, c) ^ = 1.10- 1.20 GeV, d) ^/s = 1.20- 1.38 GeV. 
The spectrum is for events, satisfying additional selection criteria: A'^^ < 7 Etot > ^.by/s. The 
shaded histogram shows estimated e^e~ — > ujtt^ background. 



initial particle. The energy dependence of the detection efficiency for the process (Q) with 
< IQ MeV is shown in Fig. |3] The decrease of the efficiency at large i/i is caused by 
degradation of resolution in kinematically fitted i^5-meson recoil mass for energies, which 
are far from KsK^ production threshold. Fig. El shows the detection efficiency dependence 
on for several energy points. 

The correction for the detection efficiency, which takes into account the difference in 
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FIG. 4: Detection efficiency for the process 

as a function of the energy of the 

as a function of energy for the events with 

photon emitted by initial particles for three 

the energy of the photon emitted by initial 

different energies: a — ^/s = 1.15 GeV, b — 

particles E.^ < lOMeV. 

= 1.27 GeV, c — = 1-38 GeV. 



detector response between data and simulation was estimated using events from the (j)- 
resonance region, where the process (0) can be separated with a neghgible background 
without constraints on x'kc;-,2tt0 d- '^^^ following selection criteria were applied: A^^ = 5, 
E^^ > 100 MeV, Z(nx5, n^g) > 130°, where E^^ is an energy of the photon not included in the 
reconstructed Ks meson, Z^nxg, n^g) is an angle between directions of the reconstructed Ks 
meson and redundant photon. For experimental and simulated events selected using these 
requirements the fractions tmc and r^xp of events satisfying standard selection criteria were 
calculated. The correction for the detection efficiency is determined as k = Texp/rMC- It is 
equal to 0.956 ± 0.015. The quoted error is statistical. 



V. DETERMINATION OF THE BORN CROSS SECTION 

The visible cross section a^s of the process under study, which is directly obtained from 
experimental data, is related to Born cross section ctq as: 

1 

CTvisiVs) = J dz- cro(v^(l - z)) ■ F{z, s) ■ e{^/s, z) (6) 
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where F(., s) ,s a probabil.tvdensity funet.ou for the initial particles to emit photon carrying 

the fraction of energy z^/s |9^], £{^/s, z) is a detection efficiency as a function of y/s and z. 
The experimental Born cross section is determined using following procedure. The measured 
visible cross section as a function of energy ams,i = Ni/ILi (here Ni is a number of selected 
events after background subtraction for the i-th energy point, ILi is an integrated luminosity 
for this point) is approximated by a function calculated using Eq. © with some model for 
the Born cross section. As a result of the approximation the parameters of this model are 
calculated together with the function R{s) = <7vis{s)/(TQ{s). Experimental values for the 
Born cross section are determined then according to the following equation: 



R{si 



(7) 



Model dependence of the result is estimated by variation of the Born cross section models. 

The Born cross section of the process e+e^ KsKl was considered in a framework of 
the VDM: 

127r 



„3/2 



E 

v=p,ll),4> 



s — my + imvTv{s) 



(8) 



The ratios of partial widths and relative phases of p, a;, and mesons were taken according 
to SU(3) model: 



Op = 0°, 6^ = 180°, 6^ = 180" 



The masses my and full widths Fy for the excitations of p, u (f> were taken as in The 
approximation is done using program FIT The following models for the Born cross 

section were considered: 



1. The process is described by four vector mesons p, u, p(1450). The relative phase 
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p(1450) 



0° 



2. The process is described by four vector mesons p, u, 

6'0(168O) = 0°, 



'1680). The relative phase 



3. The process is described by four vector mesons p, u, p(1700). The relative phase 

6'p(1700) = 0°. 
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For the models described above the following values were obtained: xl/ndf = 19.1/21, 
'X^/ndf = 18.9/21 x^/ndf = 18.7/21, respectively. All three models provide good approx- 
imation of the experimental data and can be used for the estimation of the Born cross 
section. Our final result is based on approximation by a second model. The cross section 
values are listed in Table HI Also shown in the table are the values of the detection efficiency 
and radiative corrections calculated as: 

fdz-aoiy^{l-z))-F{z,s) 
1 + 6{s) = -J- (9) 

e{s) = , ."T.^'K, (10) 
cTois) ■ {l + 6{s)) ^ ^ 

The function R{s) introduced above is defined as R{s) = s{s){l + S{s)). 



VI. SYSTEMATIC ERRORS 

The full systematic error on Born cross section of the process e~^e~ —>■ KsK^ includes 
several contributions: the integrated luminosity uncertainty, the uncertainty in the detection 
efficiency, the errors on the beam subtraction and radiative corrections estimations. 

The error on integrated luminosity. The integrated luminosity at SND detector is 
determined using QED processes e^e~ e^e~ and e^e~ 77, for which the cross sections 
are known with a precision of better than 1%. As an estimate of the systematic error we 
take the difference between the luminosities obtained using these processes, which is about 
2% almost independently of the beam energy. 

The uncertainty of the detection efficiency The analysis of the systematic uncer- 
tainty of the detection efficiency was done using events from the 0-resonance energy region. 
The cross section of the process ^ was measured using events with four or more recon- 
structed photons. The events of this class contain four photons from the Ks ~^ 27r° decay 
and additional clusters from decays or nuclear interactions of mesons or clusters from 
beam background. As it was described above the selection of the events of the process 
is based on reconstruction of Ks meson. No constrains on additional photons were applied. 
Such approach minimizes the systematic error from the inaccuracy of the simulation of the 
nuclear interaction of Kl meson in the detector. Nevertheless, since it is not possible to 
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TABLE I: Born cross section of the process e~^e~ KsKl, measured by SND detector. N^xp — 
the number of selected events, ^ Ni^^g — the number of background events, Si — detection effi- 
ciency, 1 + Si — radiative correction. The quoted errors are statistical and systematic, respectively 
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3.0 ± 1.0 


0.079 
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27.3 ± 1.8 ±0.8 


1.05 
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2.0 ±0.8 


0.078 


1.81 


15.5 ± 1.2 ±0.5 


1.06 




274 


421 


2.5 ± 1.1 


0.077 


1.92 


10.3 ±0.5 ±0.3 


1.07 




97 


96 


1.2 ±0.6 


0.076 


1.78 


7.2 ±0.8 ±0.2 


1.08 




572 


420 


6.2 ± 1.3 


0.075 


1.49 


6.4 ±0.3 ±0.2 


1.09 




94 


48 


0.9 ±0.4 


0.074 


1.30 


5.2 ±0.8 ±0.2 


1.10 




436 


158 


4.8 ±0.9 


0.073 


1.20 


4.0 ±0.3 ±0.1 


1.11 




88 


21 


1.3 ±0.4 


0.072 


1.13 


2.75 ± 0.65 ± 0.08 


1.12 - 1 


13 


420 


97 


5.6 ± 1.0 


0.071 


1.09 


2.81 ±0.30 ±0.09 


1.14-1 


15 


358 


61 


3.6 ±0.7 


0.069 


1.05 


2.22 ± 0.30 ± 0.07 


1.16 




316 


40 


2.3 ±0.5 


0.067 


1.02 


1.74 ±0.29 ±0.05 


1.18-1 


19 
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44 


4.3 ±0.8 


0.065 


1.00 


1.04 ±0.18 ±0.03 


1.20 - 1 


21 
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32 


3.7 ±0.7 


0.063 


0.99 


0.80 ± 0.16 ± 0.03 


1.22 - 1 


23 


465 


25 


3.8 ±0.8 


0.060 


0.99 


0.77 ±0.18 ±0.02 


1.24 - 1 


25 
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22 


2.2 ±0.5 


0.058 


0.98 


0.62 ±0.15 ±0.02 


1.26-1 
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27 
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20 
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1.6 ±0.4 
2.4 ±0.6 
1.0 ±0.2 


0.056 
0.054 
0.052 


0.98 
0.97 
0.97 


0.66l^:i8 ± 0-02 
QMtofj ± 0.02 
0.43l[]:i^ ± 0.01 


1.32 - 1 


33 


516 
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2.1 ±0.5 


0.050 


0.97 


o.o4j:[]:^? ± 0.01 


1.34 - 1 


35 


676 


13 


2.7 ±0.6 


0.048 


0.97 


o.33iJ]:|f ± 0.01 


1.36 




606 


11 


3.1 ±0.7 


0.047 


0.96 


0.29l^;J^ ± 0.01 


1.37-1 


38 


722 


11 


2.2 ±0.5 


0.045 


0.96 


o.28lg:if ± 0.01 
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distinguish the clusters from the i^5-meson decays from those from those produced by K^^ 
meson or beam background, these additional clusters give rise to a combinatorial background 
where lost or misreconstructed photons from Kg decay are replaced by clusters produced 
by Kl or beam background. The more such clusters, the larger number of misreconstructed 
Ks mesons. This effect leads to an uncertainty of the detection efficiency due to inaccuracy 
of simulation of the Ki-meson interactions in the detector material. This includes not only 
the inaccuracy of the total cross section of the nuclear interaction of meson, but also the 
inaccuracy of the number of clusters and energy depositions in the calorimeter counters. 

The total systematic error on detection efficiency is a sum of the uncertainty for the 
"pure" Ks meson and uncertainty of the combinatorial background increased by inaccuracy 
of the simulation of the nuclear interaction of Kl meson in the detector. 

The uncertainty of the simulation of the "pure" Ks meson was studied using events, in 
which the K^ meson is reconstructed as a single photon. The selection criteria for these 
events are described in section IIVI It was found that corresponding correction to the effi- 
ciency estimated by simulation is equal to 0.956 ± 0.015. This correction accounts for the 
differences in distributions in (3%) and photon quality parameter (2%). 

The systematic error on combinatorial background and nuclear interaction of K^ mesons 
was estimated as a difference between the detection efficiency values obtained by two ways. 
The first, standard, way is to estimate detection efficiency Smc as the ratio of the numbers 
of selected signal events to the total number of simulated events. The second way is to first 
estimate the detection efficiencies for the subsets of simulated events with fixed numbers 
of reconstructed photons and then to average obtained values according to relative weights 
observed in data. The ratio of the efficiencies obtained in this way is e\.^(j/eMc = 0.991 ± 
0.007. Corresponding systematic uncertainty is equal to 1%. 

Another source of systematic uncertainty is the difference in energy and angular reso- 
lutions for photons between data and simulation. This difference affects the resolution in 
recoil mass of the reconstructed Ks meson shown in Fig. IHl and in the slope of the depen- 
dence of the detection efficiency of the process on E^, shown in Fig. El To evaluate this 
contribution into the systematic error we varied the slope of the efficiency dependence on 
within limits corresponding to observed 2% difference in Ks recoil mass resolution in data 
and simulation on (p resonance. 

The total systematic error on detection efficiency changes from 2.1% to 2.5% with the 
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TABLE II: Systematic uncertainties of the measured Born cross section of the process e~^e 
KsKl. 



Source 


1.02 ~ 1.40 GeV 


Integrated luminosity 
Detection efficiency 
Background subtraction 
Model dependence 


2% 
2.1—2.5% 
0.1—2.9% 
0.5 — 3.0% 


Total 


2.9 ~ 5.3% 



energy in the range 1.04 - 1.4 GeV. 

Background subtraction. As it was mentioned above, the contributions of background 
processes El and IHl well agree with with estimations based on simulation. Therefore the 
extraction of the Born cross section from the experimental data the estimated by simulation 
number of background events of these processes was subtracted. The beam background was 
estimated from the x'Kg^2wO distribution of the experimental events. The systematic error 
on the background subtraction varies from 0.1% to 2.9% for the energy range from 1.04 up 
to 1.40 GeV. The statistical error on the background subtraction is included into the quoted 
statistical error of the measured cross section. 

Accuracy of the radiative corrections This systematic error includes the theoretical 
uncertainty of the radiative correction calculation, which does not exceed 0.1% jjof, and 
model dependence, related to the choice of the model describing the energy dependence 
of the cross section being measured. As an error estimate the difference between results 
obtained with the different choices of the approximation function described in section V. 
The systematic error from radiative correction uncertainty varies between 0.5% and 3.0% in 
the energy range from 1.04 to 1.40 GeV. 

The total systematic errors are listed in Table ITTl 
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FIG. 6: Born cross section of the process e^e KsKl. Dots with error bars represent results 



obtained by SND(this work), CMD-2|5|, OLYA|J] and DM1 [3] detectors. Sohd hne represents the 
cross section calculated using VDM with p{770), u;(783), (;^)(1020) taken into account. 

VII. DISCUSSION 

In Fig. inithe cross sections measured in this work and previous measurements by OLYA, 
CMD-2, and DM1 detectors are shown. The results are in good agreement with ones obtained 
by CMD-2 and OLYA in the same energy region. 

The measured cross section of the process e~^e~ KsKl significantly exceeds one pre- 
dicted by VMD model with only p(770), ^7(783), and 0(1020) resonances considered within 
SU(3) model. The curve corresponding to this prediction is shown in FiglHl The measured 
cross section significantly exceeds this VMD estimation starting from the energy about 
1.2 GeV. This excess can be described by contributions from higher mass states p', cj', and 
0', but for the determination of the parameters of these states the data in wider energy 
range and for other decay modes are needed. 

This work was partially supported by the RFFI grants 05-02-16250-a, 04-02-16181-a, 
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